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ABSTRACT 

Context. The 1(15.01 A)/I(16.78 A) emission line intensity ratio in Fe XVII has been reported to 
deviate from its theoretical value in solar and stellar X-ray spectra. This is attributed to opacity in 
the 15.01 A line, leading to a reduction in its intensity, and was interpreted in terms of a geometry 



in which the emitters and absorbers are spatially distinct. 



Aims. We study the 1(15.01 A)/I(16.78 A) intensity ratio for the active cool dwarf EV Lac, in both 
, flare and quiescent spectra. 

Methods. The observations were obtained with the Reflection Grating Spectrometer on the XMM- 

00 ! 

Newton satellite. The emission measure distribution versus temperature reconstruction technique 
is used for our analysis. 

Results. We find that the 15.01 A line exhibits a significant enhancement in intensity over the 
optically thin value. To our knowledge, this is the first time that such an enhancement has been 



detected on such a sound statistical basis. We interpret this enhancement in terms of a geometry 
in which the emitters and absorbers are not spatially distinct, and where the geometry is such 
that resonant pumping of the upper level has a greater effect on the observed line intensity than 
resonant absorption in the line-of-sight. 

Key words, radiative transfer - stars : activity - stars : coronae - stars : individual (EV Lacertae) 
- techniques : spectroscopic - Xrays : stars 

1. Introduction 

Although opacity effects have been frequently detected in the chromospheric and transition region 
emission lines of lowly ionized species in solar and stellar sources (see, for example, Christian 
et al. 2006), this is not the case for coronal transitions in more highly ionized systems. In some 
respects, this is not surprising, as the optical depth at line centre scales directly with wavelength, 
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and is inversely proportional to the square root of the plasma temperature (Mitchell & Zemansky 
1961). Coronal lines in highly charged ions will generally lie at short wavelengths (in the X-ray 
spectral region), and be formed at high values of plasma temperature. Schmelz et al. (1997) and 
Phillips et al. (1996) found some evidence of opacity in the Fe XVII 15.01 A line from an analysis 
of solar active regions. However, no evidence of coronal opacity was found in the RS CVn binaries 
Capella (Phillips et al. 2001) and II Peg (Huenmoerder et al. 2001), nor in a sample of 26 late-type 
stars (Ness et al. 2003), although Testa et al. (2004) claim the detection of opacity in the O VIII 
18.97 A emission line of II Peg. 

More recently, Matranga et al. (2005) have analysed a 50 ksec XMM-Newton observation of 
the late-type star AB Dor. They found that the observed emission line intensities were generally in 
very good agreement with theoretical values from an optically thin model. However, the FeXVll 
intensity ratio 1(15.01 A)/I(16.78 A) was significantly smaller during the flare period of observation 
than the optically thin prediction. This was attributed to opacity in the 15.01 A transition, which 
has a much larger oscillator strength than the 16.78 A feature. There was no evidence of opacity in 
the quiescent X-ray spectrum of AB Dor. 

EV Lac (GJ 873) is a very active dM4.5e flare star at a distance of 5 pc. It exhibits spot-related 
quasi-periodic variations, and has been reported to have a solar-like activity cycle (Mavridis & 
Avgoloupis 1986). EV Lac has one of the highest surface magnetic flux in the recent survey of low 
mass stars by Reiners & Basri (2007). In the present paper we study the X-ray spectrum of EV Lac, 
with emphasis on the effects of resonant scattering, via the analysis of Fe XVII line intensity ratios. 

2. Observations & Analysis 

The observations of AB Dor reported by Matranga et al. (2005) were part of a larger dataset that 
also included the star EV Lac. Observations of these stars were made during both flare and quies- 
cent periods and were obtained with the XMM-Newton Reflection Grating Spectrometers RGS 1 
and RGS2, which provide spectra with a resolution of 70-500 in the 0.2-2.5 keV range. Data were 
reduced with the XMM Science Analysis System (SAS) software v5.3.3. The RGS1 and RGS2 
spectra are divided into 3400 channels with widths ranging from 0.007 to 0.014 A, and could not 
be summed directly as their wavelength grids do not coincide, and were therefore rebinned into a 
new grid with a width of 0.02 A for further analysis. A total line spread function was defined for 
the summed spectrum, which is the sum of the individual line spread functions weighted by the 
effective area (see Scelsi et al. 2004). The analysis of the RGS spectra is based on the plasma emis- 
sion measure distribution versus temperature reconstruction technique, and was performed using 
the IDL-based software PINTofALE (Kashyap & Drake 2000). 

In our analysis, the first step was the identification of the strongest lines in each spectrum, 
using the latest version (Version 5.2) of the CHIANTI atomic database (Dere et al. 1997; Landi 
et al. 2006). The line spread function of the RGS instrument is characterized by extended wings, 
and multicomponent fitting involving the adjacent transitions was therefore necessary for most 
of the lines. In this case the total emissivity is the sum of the emissivities of the single lines, each 
multiplied by the corresponding element abundance. Lines rejected from our analysis include those 
whose inverse emissivity curves did not agree with those of other lines of the same ion. Weak lines, 
which are affected by large uncertainties and strong blending, were also excluded from the emission 
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measure reconstruction. Further details of our analysis procedures may be found in Matranga et al. 

(2005) . 

3. Results and Discussion 

The observed intensity ratios of the FeXVll 15.01 A (3d — > 2p) to 16.78A (3s — > 2p) transitions 
are consistently smaller that those calculated by collisional-radiative models. In an effort to ex- 
plain the discrepancy Laming et al. (2000) measured the ratio using the NIST EBIT and found it 
in agreement with the ratios predicted using a distorted wave collisional model. Similar experi- 
ments in the LLNL EBIT-II showed that resonant excitation can make a significant contribution 
to the 3s — > 2p line intensities (Beiersdorfer et al. 2002). Additional processes therefore have to 
be taken into account. Doron & Behar (2002) have shown that the inclusion of radiative recom- 
bination, dielectronic recombination, resonant excitation and inner-shell collisional ionization into 
a collisional-radiative model, can rectify some of these discrepancies. The Doron & Behar (2002) 
calculations and are used in our analysis. We note that the most recent calculations of Landi & Gu 

(2006) are in agreement with Doron & Behar (2002). 

In Table 1 we present the observed 1(15.01 A)/I(16.78 A) line intensity ratios for EV Lac 
during 4 different periods of observation, namely 1 quiescent and 3 flare, where the quoted errors in 
^observed correspond to 2a. The predicted ratios for an optically thin plasma from the calculations 
of Doron & Behar (2002), Rtheory < and the ratios generated from the differential emission measure 
(DEM) reconstruction, Rdem, are also given for each case. Values quoted for Rtheory are for 
2 different estimates of the plasma temperature (T e ), with Rtheory = 1-75 calculated for T e = 
7.2 x 10 6 K, determined using the measured I(Fe XVIII 16.07 A)/I(Fe XVII 16.78 A) intensity ratio 
in conjunction with the calculations of Phillips et al. (1997). By contrast, Rtheory = 1-85 is the 
theoretical ratio for T e = 10 7 K, the temperature corresponding to the peak of the Fe XVII 15.01 A 
line contribution function. 

Table 1 shows several measured 1(15.01 A)/I(16.78 A) ratios which are larger than the theoret- 
ical values, with the most pronounced example being for the quiescent phase of EV Lac, where the 
observed ratio is greater than the theoretical value by 2.6ct. (The estimated errors in the observa- 
tions are smallest for the quiescent phase because the period of observation is largest in this case). 
This result implies enhancement in the 15.01 A line intensity, as opposed to the reduction in inten- 
sity which might normally be expected for an optically thick transition. Our findings are confirmed 
in Fig. 1 , which shows the EV Lac spectra for the four periods of observation. The observed spectra 
and those predicted assuming an optically thin plasma using the DEM reconstruction are markedly 
different only for the 15.01 A line during the quiescent period. For the quiescent spectrum, the 
16.78 A line intensity is in good agreement with the optically thin model, whereas that of the 15.01 
A line is more intense than the model prediction. We would like to emphasize that the DEM anal- 
ysis of the flare and quiescent spectra included several strong lines in the 5-38A wavelength range. 
Below we give the lines used for the reconstruction of the quiescent spectrum together with their 
approximate temperatures of formation. Ne X (12.13A, logT ~ 6.8), Ne IX (13.45 A, logT ~ 6.6), 
Fe XVIII (14.21A, & 14.54A, logT - 6.9), Fe XVII (17.10A, logT - 6.7), O VIII (18.97 A, logT 
- 6.5), O VII (21.60A logT - 6.3), N VII (24.78A, logT - 6.3). The DEM distributions used in 
this work are presented in Fig. 2. 
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The statistical significance of the measurement of the enhancement depends on the degree to 
which we can safely assume that the EV Lac measurements are independent from one another 
and from the AB Dor observations (which showed no enhancement) obtained during the same 
series. We note that if they can be considered independent, the observed line ratio being larger than 
the theoretical value from the calculations of Doran and Behar (2002) by 2.6a indicates a 99.5% 
confidence level of our result being a true enhancement. If the theoretical value from the DEM 
reconstruction is used the value is reduced to 2.3a which indicates a 98.9% confidence level. We 
proceed to outline a physical scenario for such an event. 

In Matranga et al. (2005), the interpretation of the line intensity ratio 1(15.01 A)/I(16.78 A) 
being smaller than the optically thin value involved a model with a spatially separated emission 
and absorption region. To interpret this ratio being larger than the optically thin value, as found 
in the present work for the quiescent phase of EV Lac, we consider a single plasma that contains 
both emitter and absorber ions using the model developed by Kerr et al. (2004, 2005). Kerr et al. 
(2005) showed theoretically that, for spatially uniform conditions, the measured optically thick to 
optically thin line intensity ratio can be both smaller or (surprisingly) larger than the predicted 
optically thin value. Which of these occurs is dependent on two aspects of the plasma. The first is 
the overall geometry, which effectively determines the pumping of the upper state in the optically 
thick spectral line. The upper state of the transition is pumped by photons traversing the plasma 
at different angles and the characteristic length associated with this process is the mean chord 
length which is determined by the overall geometry of the plasma. Secondly, the orientation of the 
line-of-sight through the plasma determines the depletion by absorption of photons emitted by the 
upper state in the optically thick spectral line before they reach the observer. The characteristic 
length associated with this aspect is the average line of sight distance through the plasma. If the 
mean chord is longer than the average line of sight distance, then the pumping of the upper state is 
greater than the effect of depletion and the line intensity of the optically thick line is greater than 
the value assuming no optical depth. If the line of sight distance is greater than the mean chord 
length, the opposite is the case. In this way a line intensity ratio between an optically thick line 
and an optically thin line can, in principle, provide information on the plasma geometry and the 
orientation of the plasma to the observer. 

A statistical study of coronal opacity in 26 late-type stars carried out by Ness et al. (2003) using 
Fe XVII, revealed that EV Lac was the only object with significant line ratio anomalies. EV Lac has 
also been identified as one of the objects with significant optical depth in the sample of Testa et al. 
(2007). We emphasize that the previous studies were based on the intrpretation of time-averaged 
coronal spectra while in our analysis the flare and quiescent spectra have been separately analysed. 

We do not attempt here to use the observed ratio in EV Lac to determine the geometry of 
the Fe XVII emitting region in the star. This is because we recognise that, in Kerr et al. (2005), 
many simplifying assumptions were made (such as uniform plasma conditions and simple planar, 
cylindrical or spherical geometries). Instead, we note that the observed line ratio being significantly 
larger than the theoretical value is consistent with the physics elucidated by Kerr et al. (2005). 
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4. Concluding Remarks 
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To extend our interpretation requires that we develop a detailed model of radiation transfer through 
realistic stellar coronal plasmas. In addition we plan to investigate the underlying physics in two 
further ways, both of which involve measuring line ratios from plasmas of known geometry. Firstly, 
we will use the Hinode and SOHO satellites to obtain time-series spectra of optically thick and thin 
extreme-ultraviolet lines of Fe XV in solar active regions as they rotate from the solar disc centre 
to the limb. Simultaneous images in the Fe XV 284 A emission line will provide detailed informa- 
tion on the active region morphology, and in particular will allow us to construct a 3-dimensional 
picture of its geometry as it rotates from the disc centre to the limb. We will be able to search for 
and compare observed intensity enhancements in optically thick lines with the predictions of our 
models for the known plasma geometry conditions. Secondly, we believe that we can produce a 
plasma of known geometry and measurable conditions within the laboratory, for which the geo- 
metrical dependence of an appropriate line ratio can be measured. A preliminary experiment has 
been designed by Rose et al. (2006), and further work on this and similar experiments is currently 
underway. The results of both solar and laboratory studies will be reported elsewhere. Finally, we 
note that further detailed studies of EV Lac would be of great benefit in determining the level of 
confidence that could be attributed to the observations of enhancement. 
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Table 1. Observed Fe XVII line intensity ratios 1(15.01 A)/I(16.78 A), R 

observed, obtained using the 

XMM-Newton satellite for EV Lac during quiescent and flare periods. The theoretical line ratios, 
^■theory, from the calculations of Doron & Behar (2002) and the ratios determined from the DEM 
spectrum, Rdem (red line in Fig 1), are also shown. 





^observed 


^theory 


Rdem 


Quiescent 


2.50±0.50 


1.75-1.85 


1.93 


Flare 1 


2.38±1.25 


1.75-1.85 


1.77 


Flare 2 


2.13±1.00 


1.75-1.85 


1.69 


Flare 3 


1.64±0.64 


1.75-1.85 


1.89 
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Fig. 1. Observed XMM-Newton spectrum (black line), the DEM reconstructed spectrum (itself based on 
atomic data from Doron & Behar (2002)) (red line) and continuum (green line) for the star EV Lac during one 
quiescent and three flare periods. The wavelength grid of the predicted spectrum has been shifted by 0. 1 A for 
clarity. For the quiescent period only, there is a clear increase of the FeXVII 15.01 A line intensity, which is 
not observed for the 16.78 A transition. 
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Fig. 2. The differential emission measure (DEM) distribution for the quiescent state and flares of EV Lac. 
This is defined as DEM(T) = nl d ^ T ^ and is divided by 4 tt d 2 where d is the stellar distance. 
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